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Modeling strongly correlated atoms demonstrates the possibility to prepare quantum superpo- 
sitions that are robust against experimental imperfections and temperature. Such superpositions 
of vortex states are formed by adiabatic manipulation of interacting ultracold atoms confined to a 
one-dimensional ring trapping potential when stirred by a barrier. Here, we discuss the influence 
of non-ideal experimental procedures and finite temperature. Adiabaticity conditions for chang- 
ing the stirring rate reveal that superpositions of many atoms are most easily accessed in the 
strongly-interacting, Tonks-Girardeau, regime, which is also the most robust at finite temperature. 
NOON-type superpositions of weakly interacting atoms are most easily created by adiabatically 
decreasing the interaction strength by means of a Feshbach resonance. The quantum dynamics of 
small numbers of particles is simulated and the size of the superpositions is calculated based on 
their ability to make precision measurements. Experimental creation of strongly correlated and 
NOON-type superpositions with about 100 atoms seems feasible in the near future. 
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I. INTRODUCTION 

Experimental advances in controlling ultracold atoms 
have led to leaps in the ability to control and create new 
types of quantum matter [lH5(. A current challenge is 
to create larger superpositions for quantum technologies 
and for probing the foundations of quantum mechan- 
ics However, problems arise when attempting to 
make large superpositions, which improvement in experi- 
mental precision cannot overcome alone. Here, we exploit 
the quantum properties of strongly correlated atoms to 
reduce the experimental precision required. The aim is to 
provide a path for creation of the largest superpositions 
with ultra-cold atoms so far. 

Developing quantum technologies, such as precision 
measurement devices, constitutes an important goal in 
quantum physics research 0. Quantum entangled wave 
sources improve measurement accuracy without consum- 
ing more resources. Improved accuracy allows us to ap- 
ply greater scrutiny to current fundamental theories and 
has practical applications, such as gyroscopes used for 
navigation in airplanes and satellites. Current devices 
predominantly use unentangled particles, due to their ro- 
bustness against decoherence. However, some improve- 
ments have been demonstrated in gravitational-wave de- 
tectors using squeezed light [101 ] . and in atom interferom- 
etry with a Bose-Einstein condensate • If large robust 
quantum superpositions become available they will dra- 
matically improve measurement precision [l2Hl5j . At a 
fundamental level, a microscopic understanding of the 
emergence of classical behavior in large objects is still 
needed. One route to improve our understanding is by 
investigating large superpositions (lfjj . 

Currently the largest superpositions are found in su- 
perconducting loops split by one or more Josephson 
junctions, which are called flux qubits @, Hjj. They 
produce superpositions of different flux states penetrat- 
ing the loop, or, equivalently, opposite currents flowing 



around the loop. The system was initially proposed by 
Leggett [T7| and much work has been done to understand 
how such large superpositions are possible [l8rl2Cl | . Until 
recently, modeling of these types of systems relied on a 
phenomenological approach. However, by modeling an 
analogous system of ultracold atoms, previous work was 
able to show that strong interactions between particles al- 
low a large coupling between the two current states [2l[ . 
This makes the superpositions far more robust then other 
many-particle superpositions, such as NOON states (22r - 

m. 

In this paper we study several experimental factors 
that affect the creation of a binary superposition of vor- 
tex states. Initially, we study how the superposition can 
be made. According to the proposal in Ref. [2lj], exper- 
imental parameters must be evolved adiabatically after 
preparing the system in a unique ground state, such that 
the system is not excited. Two parameters are consid- 
ered here: the stirring rate of the barrier and the inter- 
action strength between the atoms. The system is most 
sensitive to the stirring rate of the barrier and we show 
how this can be increased to create the superposition. 
The fastest allowed route to formation of the superpo- 
sition is found in the strongly interacting regime. Su- 
perpositions of strongly correlated states are, thus, more 
easily accessed than NOON states, which correspond to 
weak interactions. However, it seems possible to reach 
NOON states by a different route, which is to trans- 
form the strongly-correlated superposition to a NOON 
state by adiabatically reducing the interaction strength. 
By considering the adiabaticity conditions we find that 
this second approach is much more feasible. Two other 
factors are considered: how precisely the stirring rate 
must be tuned to obtain a balanced superposition, and 
how low the temperature of the system must be kept 
for the superposition to survive. Again, we see that the 
strongly-correlated state carries the greatest potential for 
the experimental realization of a superposition. Very re- 
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FIG. 1. Ultra-cold atoms confined to a loop of circumference 
L. The position on the circumference of the loop is given by 
x — L6. The barrier is shown by the break in the ring and 
modeled by a 5-function that moves with a tangential velocity 
of v = mi ML. 



cently the non-adiabatic creation of large angular mo- 
mentum superpositions in the same system was suggested 
in Ref. US. 



II. THE SYSTEM 

A Hamiltonian describing N interacting bosonic atoms 
confined to a toroidal trapping potential of circumference 
L and stirred by a barrier can be written as 

k— — OO fci,/C2 — — oo 
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The three terms describe the kinetic energy of the atoms, 
Hk, the interaction between a barrier and the atoms, 
Hb, and the atom-atom interactions, Hi, respectively. 
Here a k creates (destroys) an atom with angular 

momentum kh. The Hamiltonian is formulated in the 
co-rotating frame of the barrier, which is assumed to be 
narrow, and is described by a 5-function with strength 
b. The tangential velocity of the barrier is v = Ml/ ML 
along the circumference of the ring, where f2 is a dimen- 
sionless rotation rate. The smallest nonzero kinetic en- 
ergy of a single particle is Eq = 2n 2 h 2 /(ML 2 ), which 
provides a natural energy unit. The radial confinement 
is assumed to be tight and has a trapping frequency of 
ujj_. The system can, therefore, be approximately de- 
scribed by a one-dimensional ring with coupling constant 
(or interaction strength) 



2h 2 a 



(a± — Ca) 



C = 1.4603. 



(2) 



w here a is the s-wave scattering length and oi = 

yj h/Mu>± characterizes the transverse confinement [28l |. 
The interaction strength can, therefore, be changed in 
two ways: the s-wave scattering length can be tuned us- 
ing a Feshbach resonance [29|, [3C| , or a± can be changed, 
which may lead to a confinement induced resonance [31] . 
This setup is schematically represented in Fig. ID 

All numerical simulations are performed for five atoms, 
N = 5, in a truncated Hilbert space with 18 angular 
momentum modes, where the interaction strength, q, is 
rescaled to account for the basis set truncation error 121, 

A general wavefunction can be written as a superposi- 
tion of terms with different total angular momentum, 



|¥>= £ C K \K) 



K=-oo 



where 



\k) = J2 (k) M 



(3) 



(4) 



is the sum over the multi-index ft = (..., no, n\, ...) 

with fixed particle number, J2k Uk = ^> an< ^ lm ~ 
plies the additional constraint ^2 k rikk — K fixing the 
total angular momentum. Using the creation operators 



a\ we can construct the permanents, 



n K 



n 



(5) 



Therefore, the probability that the state ([3]) has to- 
tal angular momentum Kh is | C^- 1 2 . In the follow- 
ing, \K), which is labelled by a single index referring 
to the total angular momentum, is a superposition of 
permanents with total angular momentum K, while 
|..., n_i, no, n%, ■■■) ...,-i,o,i,..., which is labelled by several 
indices referring to the occupation of the angular momen- 
tum modes, will describe a permanent, where nk atoms 
occupy each mode k. 

We wish to make large superpositions. One example 
is the NOON state, which is a superposition of all N 
particles being in one mode or all TV particles being in 
another mode. The NOON state can be written as 

|*)± = -^(|JV,0) 0ll ±|0,JV)o,i), (6) 

where each term in the ket represents the occupation in a 
particular mode. The terms here represent the number of 
atoms in the and h angular momentum modes. Schemes 
to create NOON states using ultra-cold atoms have been 
proposed in a double well [26|, |33( , a ring lattice of three 
sites [13, HH and a ring lattice of four sites [36[ . 

Here we consider a broader range of many-particle su- 
perpositions, which are binary superpositions of eigen- 
states of the total angular momentum operator, as op- 
posed to the two mode NOON state. This is more in 
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FIG. 2. The energy spectrum of five strongly correlated 
atoms confined to a ring trapping potential and large barrier 
as a function of Q (subsequent results use a smaller barrier 
height of b/L — 0.008i5o for increased numerical accuracy). 
At £1 = 7r there is an avoided crossing in the lowest two levels 
and this is the point where the superposition |0) + \N) is 
found. The inset shows the probability that the ground state 
has total angular momentum K = 0, -P(O), as a full line, and 
total angular momentum K = Nh, P(N), as a dashed line. 
At = 7r the probability of the two states is 0.5. 



keeping with Schrodinger's idea of a cat being in a su- 
perposition of alive and dead. In Ref. [2l| it was demon- 
strated that in systems described by Eq. (fT]) it is possible 
to make superpositions of total angular momentum of the 
form, 

|$) ± = -±=(|0)±|A0) (7) 

by applying a rotation rate of f2 = tt and having an inter- 
action strength above a certain value. Here N is the num- 
ber of atoms in the system. Therefore, the NOON state 
is just a special case of this superposition. In this pa- 
per we study the experimental control required to make 
the superposition in Eq. ([7]), but first the spectrum of 
possible states is briefly described. 



III. MULTI-PARTICLE SUPERPOSITIONS 

Non-interacting atoms - For simplicity, consider just 
the and h angular momentum modes, which is a valid 
approximation for a small barrier height. The three 
terms in the Hamiltonian ((T|) play three different, but 
crucial roles in creating the superpositions. The kinetic 
energy term, Hk, changes the energy of the two modes 
as the stirring rate of the barrier changes. At 57 — tt the 
two modes become degenerate and it is at this value that 
a superposition can form. For N non-interacting atoms 



with a finite barrier height a superposition of the form 

\4>) = - (al + al) I vac) 

1 N I AH 
= ^£^(^3^-^., (8) 

is created. This is just a product of single particle su- 
perpositions and not the state we are interested in here. 
Exact results for non-interacting atoms beyond the two 
mode approximation were derived in Ref. [21| . 

Note the importance of the barrier term in Eq. (fT}, 
Hb- This is the only part of the Hamiltonian that couples 
states with different total angular momentum. Without 
a barrier, the ground state is doubly degenerate at f2 = tt. 
With a barrier, the degeneracy is lifted leaving an energy 
gap between the ground and first excited state. 

NOON state - Equation JHJ is a superposition of 
states of the form \N — r, r)o,ij because all these states 
are degenerate for g = 0. The degeneracy is lifted 
by the interaction term, o t i(N — r, r\Hj\N — r, r)o,i = 
^ [N(N - 1) + 2r(N - r)], leaving \N, 0) and |0, N) as 
the lowest degenerate energy states. Experimentally, the 
interaction strength can be changed using a Feshbach 
resonance [29J. Again, Hb couples the |iV, 0) and \0,N) 
states and the ground and first excited state become the 
superpositions in Eq. ([5]) for a large enough interaction 
strength. 

Tonk- Girardeau superposition - For stronger interac- 
tions, the coupling to momentum modes other than 
and h becomes significant and we can no longer write the 
state as in Eq. ©. However, the ground and first excited 
states remain in superpositions of total angular momen- 
tum as described by Eq. (J7J. In the regime of strongly 
interacting bosons at low densities, and a tightly confin- 
ing wave guide, there is a one-to-one correspondence with 
spinless non-interacting fermions. The bosonic atoms are 
said to have undergone fermionisation, because they can 
no longer pass one another. This system is call a Tonks- 
Girardeau gas [37| . The energy spectrum is the same as 
that for spinless non-interacting fermions. However, al- 
though atoms cannot be in the same position, they can 
have the same momentum, so the single-particle momen- 
tum distribution is different to that of spinless fermions. 
Away from zero momentum, the distribution spreads as 
oc l/\/jk\ over an infinite range even for zero tempera- 
ture [38[. 

The energy spectrum of the Tonks-Girardeau gas is 
calculated using the single particle spectrum in the same 
way as for spinless fermions. Consider an odd num- 
ber of atoms (the calculation is more complicated for 
even numbers, but still possible 39]), where the single- 
particle energies are given by e n , £o is the lowest en- 
ergy, and the energies are in ascending order. Therefore, 
the ground state energy is then Yl n =o £ ™' ^ ne nrs * ex " 
cited state is J2n=o s u + £n, and the energy difference is 
AE = En — £n-i- Even in the Tonks-Girardeau regime 
the binary superposition of total angular momentum is 
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still created and Ref. [2l| demonstrated that this super- 
position is most robust against the loss of atoms than 
NOON states. 

Energy gap - The energy spectrum of the system pro- 
vides a useful insight into the experimental control re- 
quired to create these superpositions. Figure [2] demon- 
strates how rotation of the barrier effects the state of 
the system. At the rotation rate that produces Q = n, 
an avoided crossing is formed between the ground and 
first excited state, which is a signature of a superposi- 
tion. This can be understood by considering a simple 
two state model with states |0) and \N), which will be 
useful later. Writing the wavefunction \tp) — a\0) + /3\N) 
in matrix form, 



\1>) = 



(9) 



the effective Hamiltonian, derived from the Hamiltonian 
in Eq. is, 



H = 



A 

A E N(1 - O/tt) , 



(10) 



where A is the coupling between the two states and a 
constant energy term has been ignored. The eigenstates 
of the Hamiltonian are, 



, . _ ( c a (n) 



and jV'i) = 



C 6 (0) 
-C a (Q) 



where 



\C \ 2 = -± E o N (n-*) 

' aM 2 2VA 2 +S 2 A^ 2 (r!-7r) 2 



(11) 



(12) 



and the ener gy difference between th e two energy levels 
is hw i = y/AA 2 + E 2 N 2 (1 - Q/tt) 2 . Note that at the 
avoided crossing, ft = ir, the energy gap is just twice the 
coupling between the two states, AE = 2A. The subplot 
in Fig.[2]shows the probability of finding the system in the 
|0) and \N) states. Notice how quickly the superposition 
becomes unbalanced when Q is detuned from ir. The 
precision of tuning required to create the superposition 
is discussed in detail in Sec. IV Al 

Dependence on the interaction strength - Figure [3] 
shows how atom-atom interactions effect the energy spec- 
trum, where all energies are taken relative to the ground 
state. Notice the dip in the lowest line which signals a 
near degeneracy of the ground and first excited state. We 
find that this is the interaction strength required to make 
the NOON state. The small energy gap shows that the 
coupling between the \N, 0)o,i and |0, N)o,i states is very 
weak. The energy spectrum, together with the temper- 
ature determines the thermal state of the system, which 
is discussed in Sec. IV Bl Another factor that is strongly 
dependent on the energy gap are the conditions to reach 
adiabaticity upon changing parameters. Here we analyze 
how f2 can be evolved in time to create the superposition. 
This method turns out to only be feasible for creating su- 
perpositions in the strongly interacting regime. We study 
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FIG. 3. The energy spectrum of atoms confined to a ring 
trapping potential as a function of interaction strength g with 
b/L — 0.08i5o and Q = ir. Energies are given relative to the 
ground-state energy. 




FIG. 4. The energy level splitting at the avoided crossing, 
AE, is plotted against N for both (o) the NOON state (also 
found in Ref. [2l]|) and (b) the Tonks-Girardeau state. Part 
(a) shows analytic result for two modes (solid line) and the 
numerical result for two modes (+), where b/L = 0.008i5o and 
g/L = 47rb\/iV / 'L(N — 1). The numerical result breaks down 
for TV > 18 due to limited numerical accuracy. The analytic 
result is still valid beyond this point. Part (6) bottom figure 
shows the analytic result for b/L = O.OSEo. For large N and 
a small barrier, AE asymptotes towards AE = 26/ L. 



a more realistic method for making a NOON state that 
starts with a superposition in the strongly interacting 
regime, and then slowly reduces the interaction strength 
to that needed to make a NOON state. 

Dependence on the number of atoms - To create large 
superpositions we must also consider how the system 
changes as the number of atoms is increased. Fig- 
ure shows the energy gap between the ground and 
first excited states as a function of the number of atoms. 
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These results were obtained analytically as described in 
Ref. [2lj. The NOON state energy is calculated from a 
two mode system with the and h angular momentum 
modes, which gives, 



AE 



N 



g N - 1 L(N-l)\ 



where we require 



'N qN 
— <<: 2L <<Eo 



(13) 



(14) 



for the NOON state to be created. The energy gap, and 
therefore the coupling between the two states, rapidly 
decreases as the number of atoms is increased as seen in 
Fig. 2] (a). This makes the superposition very fragile. 

The energy gap in the Tonks- Girardeau regime is cal- 
culated using the Bose-Fermi mapping. Therefore, only 
the single particle energy spectrum is needed. At 51 = it 
this is given by the roots of 



2nh 2 a ll 
MLb 



tan(7ra A , 



(15) 



for odd fi and a M = (fx + l) /2 for even /i, where the 
energy is just e M = a^Eo [2l[. We find that the energy 
gap actually increases with atom number to 2b / L for a 
small barrier height, b/L <C NE as shown in Fig. @] (b), 
and reaches a maximum of (N — \/2)Eq/2 for a large 
impenetrable barrier, b/L > NEq. The scaling of AE 
with particle number is, thus, much more favorable in 
the Tonks-Girardeau regime than for the NOON state. 

Particle loss - Particle loss has also been shown to 
degrade multi-particle superpositions [2l|, |4(J. NOON 
states are destroyed with the loss of a single particle if the 
state of the particle is known. However, for strongly in- 
teracting atoms the effect is less detrimental. This is be- 
cause the single particle momentum distribution for the 
|0) and \N) states have a large overlap. Consequently, the 
knowledge of the angular momentum of a single atom is 
not enough to determine the constituent of the superposi- 
tion where the atom came from. The loss of a single atom 
with known angular momentum, thus, only degrades the 
superposition, but does not destroy it completely. 

Barrier shape - Another complicating factor when con- 
sidering the experimental feasibility of the setup is a finite 
width of the barrier. A barrier of finite width was consid- 
ered in Ref. [4l[ , where it was found that the energy gap 
decreases exponentially with the number of atoms. This 
can be understood simply by replacing the 5-function 
barrier by the Gaussian 



V(x) 



27r<7 



exp 



x 

2^ 



(16) 



where a defines the width of the Gaussian and b is the 
area of the potential. To first order in b and for a«L, 
the energy gap is given by 



b ( 9 er 2 9 
AE = 2- exp —2ir—N 
L \ L 2 



(17) 



which was found previously for a related system in 
Ref. [4l|. If we want the reduction in the energy gap, 
due to a finite barrier width, to be less than 1/e we re- 
quire 



a < 



V2ttN' 



(18) 



In contrast, a (5-function barrier couples all angular mo- 
mentum modes equally. 

In order to experimentally achieve efficient coupling 
the barrier width needs to be decreased with increasing 
particle number. Another solution could be to design 
a potential that can efficiently couple the relevant mo- 
mentum modes for a given number of atoms, e.g. by 
Bragg reflection off an optical ring lattice with commen- 
surate filling of atoms [25| or laser speckle. It should be 
noted that Ref. 4l| considers a different system, with 
atoms confined to a ring lattice at low filling rather than 
a continuous loop. Comparison between the results of 
Ref. [4l|, on the one hand, and those of Ref. [2l[ and this 
paper, on the other hand, suggest that the continuous 
loop is a preferable system for experimentally creating 
large superpositions. 



IV. ADIABATIC EVOLUTION 

To create the superposition of total angular momen- 
tum the experimental parameters must be evolved slowly 
enough so the system has only a small probability of be- 
ing excited out of the ground state. However, the length 
of time the system can be evolved for is limited by the 
coherence time of the atoms. Therefore, a compromise 
must be made that allows the system to be evolved suffi- 
ciently slowly so the probability of an excitation is small, 
and fast enough so the coherence of the system is not lost. 
Here we study the evolution rate necessary to achieve this 
for 57 and g. 

In general, if we evolve a quantity £ at a linear rate, 
d^/dt = constant, for a time interval [to, tf], then we can 
approximate the probability that the system is excited 
to state j by [H, HH 



Pq-jj < max 



Ce[Ca,C/] \u j) 



where hui 



oj 



Ej - E , 



(19) 



(20) 



| ^ ; and Co and C/ are the initial and final 



no = (% 

values of £, respectively. We choose the value of £ that 
maximizes the ratio aoj/woj in order to obtain an up- 
per bound on the transition probability. The probability 
that the system is not excited to state j is, therefore, ap- 
proximately 1 — e 2 . In what follows, the excitation to the 
first excited state is by far the most probable, so higher 
excited states can be ignored. 
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A. Creating the superposition 

To create the superposition, Cl must be changed from 
to 7T. For an interaction strength equal to, or greater 
than that required to make a NOON state, the superpo- 
sition in Eq. (JTJ) is created and this is what we consider 
here. The maximum probability for a transition out of 
the ground state is found at ft = tt and this is to the 
first excited state. Therefore we only need to consider 
the states |0)and \N). The effective Hamiltonian is given 
by Eq. (TTU)) . Here the coupling is calculated numerically 
using the relation A = AE/2. From Eq. (fT2"j) we find 



F(n) 



N E 

2n 2 AE' 



(21) 



Therefore, to evolve Q with only a small probability of 
exciting the system out of the ground state we need 



Po- 



< 



dfl Nh E 
~dt 2tt 2 AE 2 



< 1. 



(22) 



This shows that the system can be evolved quicker for 
smaller numbers of atoms and a larger energy gap at the 
avoided crossing. For the Tonks-Girardeau state with a 
large number of atoms and a small barrier this is 



dCl b 2 8tt 2 

~dt < E hL 2 N ' 

while for a large barrier this is 

dVt n 2 NE a 
~dt < h ' 



For the NOON state we find, 



dVt 4tt b 



2N 



(23) 



(24) 



(25) 



dt HE L 2 g 2 ^- 1 

where Eqs. (I17|) and the Stirling's approximation has 
been used, A! « \/2ttN (N / e) N . As we can see, the 
scaling for the NOON state is far worse than the Tonks- 
Girardeau state. 

B. Creating a NOON state 

We have shown that adiabatically changing O to cre- 
ate the NOON state is only possible for small numbers 
of atoms, because the near degeneracy at O = ir makes 
the evolution time much longer than an experimentally 
realizable coherence time. Another possibility is to cre- 
ate a superposition in the Tonks-Girardeau regime first 
and then adiabatically decrease the interaction strength 
to that required for the NOON state. Another advan- 
tage of this scheme is that the barrier can be removed, 
which eliminates a potential source of decoherence due 
to coupling between states with different total angular 
momentum. This must be done at a rate that does not 
change the superposition significantly, while slow enough 



not to excite the system. Now there is no longer coupling 
between states with different angular momentum. 

The properties of the ground state change markedly 
when the interaction strength is reduced from the Tonks- 
Girardeau regime to the interaction strength needed for 
the NOON state. However, we can get an idea of how 
long it will take to create the NOON state by looking at 
the most sensitive part of the evolution, which is at the 
NOON state. Because there is no coupling to states with 
different total angular momentum we will just consider 
the states with zero total angular momentum. The two 
most important states are thus |0, A, 0) _i,o,i and |1, N — 
2, l)_i.o,i, and we assume a linear ramping in time of the 
interaction strength g. 

Again, writing the wavefunction \ip) = 
a|0, N, 0)_i, o ,i + /8|1, N - 2, l)-i,o,i in matrix form 



(26) 



the effective Hamiltonian, derived from the Hamiltonian 
in Eq. flTJ|, is 



H = 



g/2L 
g/2L 2E Q + £(4N-6) 



(27) 



By calculating the eigenenergies and states of the Hamil- 
tonian we come to the solution 



F(g) 



E y/N(N-l) 



L 



(28) 



"01 



where 
huj 01 = 



2En 



-2T(4A-6) 



H2lJ 



is the energy difference of the two eigenstates of the 
Hamiltonian. 

The probability of a transition is maximised by taking 
g to zero, however, we only need to reach the interac- 
tion strength that creates the NOON state. To make the 
probability of exciting the system small we must satisfy 
the condition 



dq L 

— < ; 

dt E y/N{N - 1) 



For large N this is 



dg bVZ g 3 2 
dt HE Q l 2 



(30) 



(31) 



This scaling is far better than the scaling found in 
Eq. (|25|) . so this method represents a substantially more 
realistic method for producing large NOON states. 



V. SIZE OF THE SUPERPOSITION 

There have been several approaches to define the size 
of a superposition [l8|, [lj|, |44|. Broadly speaking these 
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fit into two categories: a measure of the entanglement 
in the system and a measure of the physical difference 
of the two parts of the superposition. Here we take a 
different approach that is based on the ability of the 
system to make a precision measurement. Superposi- 
tion states arc useful in intcrfcromctry for estimating 
parameters. Larger superpositions are better and give 
more precision. This can be quantified using the quan- 
tum Fisher information (QFI), Fq, which describes the 
maximum amount of information that can be extracted 
about a certain quantity and is the standard parameter 
used to determine the usefulness of a system for mak- 
ing a precision measurement [45l - l47l |. This is related to 
the uncertainty of the measurement via the Cramer- Rao 
bound, 8<j> > 1/ \/Fq [4^]. It has been shown that large 



superpositions can be used to perform precision measure- 
ments beyond the classical limit EH • Here the QFI 
is calculated to determine how the superposition is de- 
graded under non-ideal experimental conditions. The in- 
fluence of particle loss on the QFI was already considered 
in Ref. 

Consider a superposition of two states with different 
total angular momentum, 



|*> = ^(|#i) 



(32) 



A sudden change in the rotation rate of the barrier will 
change the energy of the two states and they will acquire 
different phases, 



f 



\Ki) + e 



-1K1 



l#a» (33) 



where a global phase factor due to the interaction and 
barrier terms have been ignored. The phase <f> — 
httt/ML 2 is the rotation angle of the barrier acquired 
during a time t by time evolution with the Hamilto- 
nian (fTJ). From the 4> dependence of the wavefunction, 
the QFI can be calculated, which measures how well the 
state can determine (p. The QFI is given by 



Fq = Tr[p(^>)A 2 



(34) 



where p((j>) is the density matrix of the system, and A is 
the symmetric logarithmic negativity, defined as 



(35) 



In the eigenbasis of p((j>) this is (A)^ = 2[p' ((j>)]ij / (Xi + 
Xj), where X L j are the eigenvalues of p((j>) and p'{4>) — 
dp(<f>)/d(f>. If Aj + Xj ; = then (A)^ = 0. For a pure 
state, \V ((/))), the QFI is 



(36) 



where \t>'(<p)) = d\V(<j>))/d<f> g|. 

For the superposition of total angular momentum 
states l|32p. the quantum Fisher information is (Ki — 



K2) 2 ■ For the case of non-interacting atoms described 
in Eq. © we do not have a binary superposition, how- 
ever we can still calculate the QFI using Eq. (|36j) and 
obtain Fq = N. This is called the shot noise limit and is 
the maximum QFI achievable by unentangled particles or 
independent measurements. A significant increase in the 
QFI is achieved when using NOON states, such as the 
one given in Eq. ©. The QFI is now Fq = N 2 . This is 
called the Heisenberg limit and is the maximum precision 
predicted by quantum mechanics [49l - l52l |. Furthermore, 
all states that can be written in the form of Eq. ([7]) also 
have this QFI. This was first demonstrated in Ref. [40| . 



A. Potential decoherence 

A major difficulty in creating large superpositions is 
their extreme sensitivity to experimental inaccuracies. 
This was first investigated for NOON states created in 
a ring lattice [25| and subsequently for NOON states cre- 
ated in a double well, where this effect was aptly named 
potential decoherence [26| . In this section, the investiga- 
tion is not limited to NOON states, but instead looks at 
the range of superpositions from uncorrelated to strongly 
interacting atoms. The investigation is restricted to im- 
perfections in the rotation rate of the barrier, because 
the system is most sensitive to this parameter. 

The same effective Hamiltonian described in Sec. IIV Al 
and given by Eq. (|10[) can be used here. Consider inter- 
action strengths strong enough to create a NOON state 
and stronger. This ensures that the ground state is a 
binary superposition of total angular momentum when 
il = it, so only the states |0) and \N) are relevant. To 
measure the effect of detuning of f2 away from ir we study 
how the QFI changes. 

Again, the coupling term is given by half the energy 
gap at Q — 7T, and can be obtained numerically. The 
eigenvalues and states of the Hamiltonian are given in 
Sec. nvXl thus 



F 



Q 



(r/2) 2 N 2 



(r/2) 2 + (n-n) 2 



(37) 



This is a function of f2 with the shape of a Lorentzian 
distribution 15311 . width 



r = 



nAE 



(38) 



and a maximum amplitude of N 2 . Equation Q36[) is sim- 
plified to F Q = AY^k [K 2 \C K \ 2 {\ - \C K \ 2 )} for the two 
states considered here. We have calculated the QFI nu- 
merically and plotted it against the analytic result in 
Fig. [5] The numerical and analytic results give good 
agreement except for small interactions. This is because 
the two state approximation fails in this regime. For 
large interactions, g, the value of the QFI is A^ 2 . How- 
ever this drops to for small g showing that the system 
is uncorrelated (see Fig.[S]). 
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0.015 
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0.005 
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FIG. 5. The inset shows the QFI as a function of Q in the 
Tonks- Girardeau regime. Shown in the inset is the full width 
half maximum, T. In the main figure, F is plotted on the left 
axis against the interaction strength, g, for for b/L = 0.08-Bo. 
The solid line shows the numerical result and the dashed line 
the semi-analytic result given by Eq. (|38[) . The right axis 
shows the QFI for different interaction strengths when = tv 
(dotted line). 



r must be large to create states that are less sensitive 
to changes in fi. This is realised when the energy gap is 
large and the number of atoms is small. AE is approxi- 
mately constant as a function of atom number for large 
numbers of atoms in the Tonks-Girardeau case, while it 
decreases rapidly for the NOON state (see Fig. @|. The 
gradients of the energy levels as a function of f2 are pro- 
portional to N, which explains the number dependence 
in the denominator of Eq. ()38[) . This will present a limit- 
ing factor for experimentally realizing this state for large 
numbers of atoms. 



B. Non-zero temperature 

The final factor considered here is how temperature 
affects the state. For a canonical ensemble in thermal 
equilibrium at temperature T, with Hamiltonian H, the 
density matrix is [53| . 



P(T) 



exp(- 



H 

k B T 



Tr 



exp( 



(39) 



where k B is the Boltzman constant. Again, the system 
can be understood by considering the effective Hamilto- 
nian given by Eq. (fit))) . It can be shown, using Eq. (fM|). 
that the QFI is 



F Q = N 2 tanh 2 



AE 

2k B T 



(40) 




Figure [B] shows the QFI as a function of temperature 
for different interaction strengths. The thin (black) lines 



k B T/E 



FIG. 6. (Color online) QFI as a function of temperature 
and interaction strength. Each line represents a particular 
interaction strength. The thick blue line show the analytic 
result for a Tonks-Girardeau superposition of Eq. (|40p and 
the red crosses show the QFI at T = from Fig. [5] The 
barrier height is b/L = 0.08-Eo- 



show the numerically calculated QFI. This was done us- 
ing Eq. and the lowest 20 energy levels were used 
to create the thermal density matrix. The density ma- 
trix was also created with larger numbers of energy levels 
to check the convergence, which showed that there was 
no significant change to the QFI. As expected, the QFI 
decreases as the temperature increases, because higher 
excited states become populated. However, the decrease 
in Fq is less rapid as the interaction strength is increased. 
The thick (blue) line in Fig. [B] shows the analytic result 
for the Tonks-Girardeau case, which agrees well with the 
numerical results. For the NOON state the decrease is 
so rapid that Fq w on this scale. The red crosses in 
Fig. |5] show the QFI at zero temperature. For interac- 
tions less than that required to make a NOON state, the 
energy gap between the ground and first excited state 
increases, so the system becomes less affected by tem- 
perature. However, the QFI at zero temperature is also 
reduced due to the loss of the binary superposition. 

For a superposition in the Tonks-Girardeau regime, 
as the number of atoms is increased, the superposition 
becomes more robust to temperature, however, this is 
due to the increase in the energy level splitting between 
the ground and first excited state. As shown in Fig. |4l 
the splitting becomes almost constant for more than ten 
atoms. This suggests the robustness to temperature will 
reach a constant value for all numbers of atoms and a 
small barrier, however this robustness can be improved 
by increasing the height of the barrier. 



VI. EXPERIMENTAL REALIZATION 

In this section we consider the experimental realization 
of 100 7 Li atoms confined to a tight ring trapping poten- 
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tial of radius 50/im. This gives a mean atom spacing of 
3.1/im. Here we choose a modest transverse oscillator 
frequency lu± — 9 kHz, which gives a transverse radial 
confinement of aj_ = 1/im. The superposition is most 
readily created for a large barrier height, however, this 
strongly couples states other than |0) and \N). As a com- 
promise we choose half the limiting value for the energy 
gap in the Tonks- Girardeau regime, AE w 25Eq m 45fi. 
Hz. 

The superposition is created by adiabatically increas- 
ing the angular velocity of the stirring barrier to 0.29 x 2tt 
Hz. This should be done with a barrier with a width less 
than 0.5/mi, according to Eq. <jTSj) , so that the energy 
gap at the avoided crossing is not significantly effected 
by a finite width barrier. From Eqs. ([37| and ([38]) we 
see the accuracy of the angular velocity must be within 
0.01 x 2tt Hz to achieve a well balanced superposition. 
From Eq. (|24|) we also see that the rate at which the ro- 
tation can be increased, with only a small probability of 
excitation, must be much less than 26 x 2it Hz s _1 . If 
we take the extreme case that the atoms are in thermal 
equilibrium with the environment, then, from Eq. (|40|) we 
would have to reduce the temperature of the condensate 
to less than 0.2nK. This is due to the small energy gap 
at the crossing. However, if we assume the population 
in the ground state remains the same as when the con- 
densate is created for a stationary barrier, then higher 
temperatures should still produce the superposition. Fi- 
nally, to create the NOON state, Eq. (|3"Tj) shows that the 
interaction strength, g, must be reduced at a rate much 
less than 1.7 x 10 -39 kg m 3 s -3 to avoid excitations, which 
is equivalent to reducing the scattering length at a rate 
much less than 0.0044 A(ms) -1 . We should note that 
this is the rate required at the value of g « LEq/(2N) 



where the NOON state is reached and excitations are 
most likely. Much faster rates will be possible for larger 
interaction strengths. 



VII. CONCLUSION 

We have studied a spectrum of multi-particle super- 
positions of total angular momentum to understand how 
to make them large and experimentally accessible. Com- 
pared to NOON states, the results presented here show 
that superpositions using strongly correlated atoms are 
far more robust to potential decoherence, finite temper- 
ature and allow a faster evolution of experimental pa- 
rameters without a significant probability of excitations. 
Although it looks unlikely that a large NOON state can 
be created by changing the stirring rate of the barrier, 
a more realistic approach would be to create a super- 
position in the strongly correlated regime first and then 
slowly ramp down the interaction strength. This could 
be achieved via a Feshbach resonance. Experimentally 
realising this scheme will provide a means to probe the 
boundary between quantum and classical mechanics, and 
improve precision measurements beyond the shot noise 
limit. 
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